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The vibrational spectra of glass formers follow different laws with respect to crystals. A rationale 
for their anomalous behaviour is provided by the euclidean random matrix theory. Experiments on 
glass formers at different densities might be a suitable route to test the prediction of the theory. 
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I. INTRODUCTION 

The high frequency (Thz) dynamics of the density fluctuations in highly viscous supercooled systems is presumably 
due to purely vibrational modes. However, with respect to their crystalline counterpart these systems show a number 
of puzzling feautures whose theoretical interpration is still controversial (see e.g. The purpose of this work is to 
focus on two of them describing a theoretical approach which is able both to explain the physical processes beyond 
such features and to make rigourous predictions. 

Brillouin-Peak X-rays and Neutron inelastic scattering experiments probe excitations of momentum q comparable 
to the characteristic momentum qq, the position of the first diffraction peak in the static structure factor. In 
other words the wavelength is of the same order of magnitude of the inter-particle distance and an hydrodynamic 
approximation does not make much sense. Nevertheless the position w and the width T of the Brillouin Peak 
are reported to show hydrodynamic-\ike features: 

Wpeafc oc q, T peak ccq 2 (1) 

The agreement regarding the exponent 2 in the scaling law of T between hydrodynamics and experiments is 
likely fortuitous. In fact the magnitude of T has a very slight dependence on temperature (if any) while it should 
be strongly temperature dependent in case the decay of hydrodynamic modes was induced by the viscosity Q . 

Boson Peak Raman experiments and neutron scattering experiments at large q (see below) may be useful to deter- 
mine the vibrational density of states (VDOS) g(w). Glass formers are quite peculiar with respect to crystals 
since in the frequency region where the main hypothesis of Debye holds (linearity of the dispersion relation) one 
finds that the VDOS has a different behaviour from the Debye law: 

g(u)/gDebye(oj) » 1. (2) 

This excess of states over the Debye VDOS is termed Boson Peak (BP). 
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II. THEORY 

Within the harmonic framework one assumes that the VDOS is simply given by the eigenvalues of the Hessian 
Matrix /Cij ;AI „. At low enough frequency this approximation breaks down since the relaxational modes start playing a 
relevant role. One could roughly claim that the harmonic approach holds whenever to is larger than the inverse of the 
fastest relaxation mode time-scale. Since this quantity in all the realistic cases is unkonwn, one must content to assess 
empirically the range of reliability of the harmonic approach by means of the comparison between the experimental 
spectra and the numeric spectra [10| . Aiming to modelize the vibrations in glass- formers one has to take into account 
the role played by the disorder. As a matter of fact the positions around wich the atoms vibrate harmonically do 
not form in general a regular lattice. The disorder of the positions {xi} of the center of oscillation is described by a 
suitable distribution of probability Different choices for P describe different physical situations. Furthermore 

this topological disorder induces a broad distribution of probability of the spring constant K,. The theoretical effort 
is greatly simplified if the macroscopic quantities are assumed to be self-averaging, i.e if their value remains the 
same when different realizations of the disorder are taken || . This allows to replace the computation of (say) the 
VDOS for a given realization of the disorder (involving the knowledge of an infinite number of positions Xi) with 
the computation of the averaged VDOS (which involves the knowledge of P). Besides the VDOS another interesting 
quantity is the Vibrational Dynamic Structure Factor (VDSF) S{q,uj), which shows the Brillouin Peak. The two 
quantities are written as: 
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where e n are the eigenvectors and lo 2 the eigenvalues of the Hessian matrix /C. The above expression of the VDSF 
arises when only 1-phonon process are considered. In this work we address the computation of such quantities for 
topologically disordered models, which we claim to be the most appropriate to describe glass formers. 
In order to compute the spectral properties of the matrix /C 

A'.. ,:.;„. = S ij X! 9 ^ V (\^k\) - 9 M i/V(|fjj |) fij = Xi - x 3 . (4) 

fc 

we introduce a complex valued tensor G^ called resolvent, which can be split in its longitudinal and transversal parts: 
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z is a complex number related to the value A of the eigenvalues of the Hessian matrix K, by z = A + ie. 

The resolvent is related to the VDOS and the VDSF. However G M „ is defined in the space of eigenvalues which 
nothing prevents from being negative. On the other hand the VDSF and the VDOS are functions of the frequency, 
and the frequencies are related to the eigenvalues only in the positive eigenvalues region (A = lu 2 ). Then the following 
relations hold only in that region: 

2kTq 2 to 2 
S(q,u) = lm.G L (q,uj 2 +i0+), g{w) = lim — -^S(q, lu) (6) 



However we will show that the existence of the resolvent at negative A leads to important conseguences. Note that 
only isotropic modes survives in the VDOS because G llv (q = oo, z) cx 5^ v . Thus the polarization is approximately 
defined only for q ~ 0. Ihis is very different from lattice models, where g(co) cx J dqS(q, u>). 

The resolvent can be computed in a non perturbative way if the n-th moment of P is assumed to be factorized 
(superposition approximation): g^ n \xi . . .x n ) — g(r 12 ) ... g(r n -in) where g(r) is the pair probability distribution. 
With this approximation the resolvent is given by [3j: 
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where p is the density, V^„(fc, q) = p \J^u{k)) ~ f^v{q ~ k)\ and f^{q) = J-[g(r)dfj_ u v(r) (v(r) is the pair interaction). 

The quantity £ is called self-energy and describes the loss of energy of phonons due to the disorder. It vanishes in 
the limit of infinite density, where phonons propagate without dissipation. In this limit one has: 



S(q,u) = 6(u>-u>° eak (q)) ^(q) = yj p (f(Q)-h(q)) 



(8) 



This does not imply that the disorder vanishes when p — > oo. Rather his effects on the propagations of phonons are 
averaged out by the infinite number of particles within each wavelength. On the other hand at finite p the self-energy 
starts playing a relevant role, modifying the position and the width of the peak. 



III. EVOLUTION OF THE BOSON PEAK WITH DENSITY 



Solving the integral equation J2J in the limit of q — > oo one obtains the density of eigenvalues from which the VDOS 
g(uj) is obtained. Due to the complexity of the equation, the result for a given choice of fuu(o) can be computed 
only numerically. However it is possible to deduce some important mathematical properties of the solution which 
are not dependent on / (hence, on the model). Keeping fixed the other thermodynamic parameter, it can be shown 
that at low enough densities the Hessian matrix has both negative and positive eigenvalues, while at high densities 
all the eigenvalues are positive. A critical density p c separates the two regions. In the low density regime © implies 
that gipS) cx u> when uj ~ while at very high densities one has the standard Debye behaviour g(u>) cx lo 2 . The 
change of exponent in g(uS) is due to the occurence of a sort of phase transition. The parameter which induces 
such phase transition is not necessarily the density, a changement in any of the other thermodynamic quantities may 
lead the system from the region with only positive eigenvalues of K, to the region where also the negative ones are 
allowed. For example it has been shown 6] that the value of the potential energy of the stationary points (minima 
or saddles) around which a glass system remains for a long time (termed generalized inherent structures) are one of 
those quantities. For the experiments however the stationary points are not a suitable observable, while it is surely 
conceivable perform measurements of spectra at different densities putting the system under a very high pressure. 
In the following we will identify the region without negative eigenvalues with the glass forming liquid. Let us see in 
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FIG. 1: Left: The VDOS in log-log representation compared with the u> 2 and lj 3 ^ 2 laws predicted by the theory. Right: The 
position of the Brillouin Peak as a function of the momentum q. 

detail the behaviour of the VDOS in this phase. In the region of frequencies where the dispersion relation is still 
linear one finds the following universal form: 

! x 2 -t x <C 1 
(9) 
const. i> 1 

If we define a distance A from the critical point as A = (p — p c ) we see that at a charateristic frequency lobp it occurs 
a crossover from the Debye behviour g(u>) ~ u> 2 to a different law g(u>) ~ w 7 . One obtains furthermore that: 

LOBP ~ A", g(uJ B p)/gDebye(uBp) ~ A~ ? ' (10) 

where the universal exponents are: 
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The main conseguence of the a critical point is the modification of the Debye law at to bp- This is higly remindful 
of the Boson Peak found in experiments. This suggests to identify to bp with the position of the Boson Peak and 
g{u>Bp) I gDebyei^Bp) with its height. Since QjJJ does not hold at low u> the predicted singularity at p c is not approached 
in real systems. However the proximity of such singularity is enough to break the Debye law. 

We can test our theoretical prediction by performing numeric simulations of a simple glass former. We performed 
MonteCarlo simulations at different temperatures and densities of a simple glass former (see pfl of the details). The 
VDOS's have been computed by diagonalizing at least 50 different realizations of the Hessian matrix. 

In fig we show the VDOS for T = 0.897T mc (T mc is the Mode Coupling temperature []]). At uj B p ~ 1 (where 
the dispersion relation is still linear) there is the crossover in the VDOS from the Debye law to the w 3 / 2 behaviour. 
This is the signature of the Boson Peak, as predicted by the theory. In fig @ we see the large variations of the VDOS 
for different densities at T — 4.420T mc (left). The position of the Boson Peak moves linearly with the density while 
its the growth of its eight seems to be compatible with a power law whose exponent is 1/2 (right). Thus even those 
theoretical predictions are in reasonable agreement with the numeric data. 
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FIG. 2: Left: VDOS divided by lo 2 at different p. Right: The position and the heigth of the BP as a function of the density 



IV. THE COMPLEX BEHAVIOUR OF THE BRILLOUIN PEAK 



At finite densities the propagation of phonons is affected by the loss of energy due to the disorder which make 
the time-life of phonons finite. Actually it is not entirely correct to describe the effect of disorder on the phonons 
in topologically disordered systems in terms of phonon-disorder interaction. In fact as we said above in the infinite 
density limit the phonons propagate freely even if the positions of the particles remain disordered. Actually, when 
the disorder has a topological nature it is not possible to find a limit where the disorder vanishes, even if his effects 
may disappear. This is at odds with lattice systems, where phonons propagate freely only when the disorder is not 
present. The position and the width of the Brillouin peak are affected by the self-energy in the following way: 



Lpeak{q) ~ 



(12) 
(13) 



Aiming to obtain model independent results we limit to the asymptotic region where q is small with respecto to qo- 
The general result is: 



u P eak(q) ~ cq 

t-i I \ A di^peak) 2 , r> 4 

Tpeakiq) ~ A —q +Bq 



(14) 
(15) 



(c is the speed of sound and A, B two model dependent constants). While the eq. (|14|) signals trivially the existence of 
propagating phonons, eq. i|15|) reveals a quite rich behaviour of the width of the Brillouin peak. When the frequency 
of the Brillouin Peak lie below ujg p the leading term of T pea k is oc q 2 while at larger frequencies there is a crossover to 
r pea fe oc g 3 / 2 . Then ujbp signals not only the frequency with the maximum excess of states with respect to the Debye 
law but even a crossover in the power law of the width of the Brillouin Peak. Moreover at large enough momenta 
there exists another crossover frequency where the q A term becomes dominant. 

We have seen in the previous section that for very large densities the BP moves a very large frequencies hence it is 
reasonable to expect that in this situation one should observe only the crossover at from the r oc q 2 to Y oc q 4 . Not 
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only lobp changes with the density (see eq. 1)11)0 h even the characteristic frequency does, though with a different 
behaviour. At high densities A cx 1/p 2 , B cx 1/p hence uu* cx at p >> 1 |3J. 

V. CONCLUSIONS 

The euclidean random matrix theory provides a coherent description of the anomalous features found in the vibra- 
tional spectra of glass formers. It provides furthermore well defined quantitative laws which should be compared with 
the numeric and experimental findings. In this paper we presented mainly the laws describing the evolution of the 
Boson Peak and of the Brillouin Peak when the density of the system is changed. 

How many chances do the experiments have to investigate such laws? The main restraint we see is given by 
the limited frequency windows where the harmonic modes describe entirely the dynamics. In fact already in the 
Ghz region the relaxational modes (accounted for example by the Mode Coupling theory) could become important. 
Optimistically we might argue that the harmonic approach holds at most over three decades (say for frequencies 
S [0.01 — \{)]Thz). This raises some problems in the verification both of eq. (|10H and of eq. 1)15(1. In the former 
case one cannot get very close to the critical point underlying the Boson Peak since the low frequency region is not 
described by harmonic effects. Still, the numeric results in fig. J2J) suggest that the singularity might be discerned 
and studied quantitatively even in experiments. In the latter case the possible presence of two crossovers is hard to 
detect having only three decades of frequencies available. Not surprisingly different experimental groups interpret the 
results in this region in different ways [2J, |4| . In that kind of meausrement it might be useful to remind that playing 
with density is it possible to change the values oi lobp and exalting a particular behaviour with respect to the 
others. 
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